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The collinear QCD structure of the electron is studied within the Standard 
Model. The electron structure function is defined and calculated in lead- 
ing logarithmic approximation. It shows important contribution from the 
interference of the intermediate electroweak bosons. The problem of mo- 
^f-) ■ mentum scales is extensively discussed. The master equations for the QCD 

\Q ' parton densities inside the electron are constructed and solved numerically 

f^) , in the asymptotic region. Significant corrections to the naive evolution 

• procedure are found. Phenomenological applications at present and future 

\ momentum scales are discussed. 
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1. Introduction 



It is the nature of the Quantum Chromodynamics (QCD) and its run- 
ning coupling constant which cause point-like objects of the electroweak 
theory to develop quark and gluon structure. This QCD component was 
first calculated Q in the case of the photon. In the forward direction, when 
(nearly) massless quarks are emitted collinearly from the photon, a consec- 
utive QCD cascade developes practically 'at no cost'. The presence of such 
a 'resolved' photon has been observed in many experiments ||j allowing now 
for a precise study of its dependence on the momentum fraction x and the 
momentum scale P 2 . Analogous calculation extended to the case of the weak 
bosons W and Z has also been performed ||. It shows several differences as 
compared to the photon case. Due to the nature of weak couplings the QCD 
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parton densities depend strongly on the quark flavour and in most cases the 
polarized densities turn out to be nonzero. The longitudinal bosons show 
much weaker structure and in the leading logarithmic approximation can be 
neglected. Phenomenological applications of the 'resolved' W and Z bosons 
require momentum scales much higher than their masses. This is because 
only in this region the approximations used are valid, moreover, only there 
the photon contribution does not dominate entirely the weak bosons' effects. 

In order to study the boson structure functions we need a source of real 
electroweak bosons. But even in the case of the photon, 'real' means in most 
cases only 'nearly on-shelk. The situation is even more approximate with 
the weak bosons. Here 'nearly real' means in practice that their momentum 
squared is negative and close to zero. This is because the best known source 
of high energy electroweak bosons is a leptonic beam. The procedure applied 
usually is to replace the incoming lepton with the spectrum of equivalent 
bosons []|, [|, Fg , and to convolute it with the parton distributions, F^, 
inside the boson, which results in the following QCD parton h distribution 
inside the electron: 

Ff(z,Q 2 ,P 2 ) = ]T [dxdy5(z-xy)F*(x,P 2 )F e B (y,Q 2 ), (1.1) 

B=7,W,Z J 

where P 2 is the hard process scale, Q 2 is the maximum allowed virtuality of 
the boson (usually taken to be proportional to P 2 ) and z is the momentum 
fraction of the parton h with respect to the electron (detailed definitions 
follow) . 

The procedure presented above has, however, several unclear points. It 
does not answer the questions: what is the relation between the two entering 
momenta scales Q 2 and P 2 , how far off-shell can the electroweak bosons be, 
are there any interference effects between the intermediate bosons (is the 
sum actually diagonal in B and/or in the boson polarizations). One may 
also ask the question: when is the convolution Eq.(|l.l|) justified and what 
are the corrections? The best way to clarify these doubts is to address a 
direct question: what is the QCD content of the incoming electron? In this 
paper we demonstrate how the concept of the electron structure function 
allows for a detailed study of the QCD cascade originating from an electron. 

The paper is organized as follows. In Section 2 we discuss the general 
structure of the QCD master equations in the system of leptons and elec- 
troweak bosons. The leading order calculation of splitting functions of an 
electron into a quark (antiquark) is demonstrated in Section 3. The problem 
of momentum scales and interference effects are discussed there in detail. 
The constructed master equations for the parton densities inside the electron 
are solved in the asymptotic region in Section 4. We demonstrate there the 
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flavour and spin dependence of the asymptotic distributions. Comments on 
phenomenological applications, summary and conclusions close the paper. 



2. Master equations 

In order to study the QCD content of a lepton £ we look for partons 
produced by a highly virtual gluon G* scattering off £: 

G* + £—>£' + h + anything . (2.1) 

where h = quark q, antiquark q or gluon G and £' is the outgoing lepton. 
We have chosen the gluon G* as a probe because it does not couple directly 
to the lepton. Our process, Eq. (|2.l| ), is a part of many physical reactions, 
the far off-shell gluon may originate from an incoming hadron or another 
lepton. The virtualness P 2 of this gluon sets the scale for the process. 

Before going to a detailed calculation we present a heuristic derivation of 
the master equations for parton densities involved in the process fl2,ip . We 
point out which assumptions usually made in this scheme are not fullfilled 
in a general case. The standard master equation for any 'parton' A density 
inside £, F^, reads 0] 



d 
~dt 



F\{z,t)=Y J I dxdy6(z-xy)V^(x,t)F^y,t) (2.2) 
B 

or, using a shorthand notation for the convolution 

4^(*) = E^(*)® i! 3S(*)- (2-3) 
at B 

In general t = \og(P 2 /Pq ), where Pq is some momentum scale. We discuss 
the choice of Pq in the next section. The 'splitting function' V®(x,t) is 
the probability density per unit t to find A carrying fraction x of i3's mo- 
mentum in a collinear decay of B. The above master equations are quite 
general within the approximation that the probability to produce a parton 
in a fusion process is negligible. It is assumed here that this probabilistic 
approach holds for a physical process i.e. the contributions from different 
'intermediate' partons B do not mix. As we will show in the next section, 
this is not the case for the electroweak sector. 



Since there is no direct lepton-gluon coupling the process (2J-) must be 



mediated by electroweak bosons and we have to consider following types of 



4 



'partons' inside lepton I: q, q, G, £, 7, W 1 * 1 , Z°. Assuming for the moment 
that the probabilistic approach holds, Eq. fl2.3|) takes the form 

£' S=7,W,Z 

+ E ( 2 - 4 ) 

ft=q,q,G 

Within the approximation of the lowest order in the electromagnetic 
coupling a and leading-log order in the strong coupling a s (t) 

F$,(x,t) = Sf 5(1 -x) + 0(a 2 ) , (2.5) 
VUx,t) = ^(x), (2.6) 

V?(z,t) = ^if(x), (2.7) 

= (2.8) 

and all other P§ = 0. 

Inserting above relations into Eq.(2.4) we arrive at 

S F » (t) = S E + ^ E (2.9b) 

S=7,W,Z h'=q,q,G 

The solution to Eq.(2.9a) is known as the density of 'equivalent bosons' 
in the lepton £ §. The splitting functions vanish for h = G. P^ and 
PjP have been calculated for B = 7 in Ref. and for B = W, Z in Ref. 

|J. P^' are the standard Altarelli-Parisi splitting functions Q. Inserting 
equivalent boson density Pg into Eq.(2.9b) one can solve it for Ffc. Eqs.(2.9) 
suggest also that there is one common momentum scale governing the whole 
system. This is not true in general as will be demonstrated below. 

As already stated the above reasoning neglects interference effects, here 
possible between 7 and Z bosons. Because of different masses and different 
couplings of these bosons to leptons this interference is non-trivial, i.e. can- 
not be removed by any kind of 'diagonalization' procedurecl. Thus, instead 
of Pg, one should rather use a density matrix, F^ B , of equivalent bosons 

The Standard Model gauge group is not simple as opposed to QCD where all colour 
interferences are taken into account by a group-theoretical weight factor independent 
of the couplings and kinematical variables. 
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inside the lepton I. In the following section we perform an explicit calcula- 
tion which shows how this density matrix arises and how to generalize the 
first term of Eq.(2.9b) in order to take the 7-Z interference into account. 

The use of the density matrix implies — at the first sight — that the 
probabilistic approach expressed in terms of the master equations breaks 
down. If, however, instead of looking at the contributions from different 
electroweak bosons we introduce the notion of the hadronic content of the 
lepton, the probabilistic interpretation is recovered and the master equations 
read 

/i'=q,q,G 

where Pl(t) is 'lepton to hadron splitting function'. In the next Section we 
define and calculate this function. 



3. QCD current of electron 

The hadronic structure of a lepton can be analyzed by looking at any 
process where the lepton couples to the QCD current. In the following we 
will calculate the matrix elements squared of this current in the case of 
£ = e~. 

In the lowest order in the electromagnetic and strong coupling constants 
(a and a s ) the current matrix element squared for polarized particles reads: 

J^ihP) = ^J dT v dT k dT k ,(2it) A 5 A {k + k' + l l -p-l) 

x E E( e Al^(o)l4'q.q^)(^q.q-,l^(o)|e^), (3.1) 

e'=e,v X'=± 

where 

d 4 k 

dT k = j—^2<K5{k 2 ) (3.2) 

and indeces A,A',ry denote helicities of corresponding particles. We do not 
sum over quark flavours and helicities because we are interested in the den- 
sities of polarized partons in the electron. In other words we choose a single 
flavour of particular helicity in the final state. 

The Feynman graphs contributing to our process are shown in Figure]!]. 
The incoming electron e _ carries 4-momentum I and the off-shell gluon G* 
of 4-momentum p with large P 2 = —p 2 , supplies the QCD current. In the 
final state we have massless quark q and antiquark q of 4-momenta k and k! 



6 



and lepton I' (electron or neutrino) of 4-momentum V '. The lepton interacts 
with the quark exchanging an electroweak boson B = 7, Z, W which carries 
4-momentum q (Q 2 = —q 2 ). For massless quarks J${l,p) is known to have 
only 3 independent components for any given polarizations A, rj. They can 
be conveniently paramatrized by the helicities of the incoming gluon. Let 
us chose a reference frame where the gluon momentum is parallel to the 
z-axis, p^ = (po,0,0,p z ). Its polarization vectors e^(p) read 

e£ = i=(0,l,±i,0), (3.3a) 

e oiP) = ^(P*,0,0,po), (3.3b) 

where P = \/\p 2 \- 

The three helicity components of J 

JF(!,P) = eftfo) J#(!,p) <fo(p) (3-4) 

are proportional to the cross-section for the (virtual) process 

e^ + G* 4 + q„ + q_ T; . (3.5) 

We will proceed with the calculation of in the reference frame where 
the B boson momentum q = —p, i.e. in the center-of-mass of B and G*, 
which is also the rest frame of the qq pair. We choose z-axis in the direction 
of q, the transverse directions are measured then in the xy plane. In this 
frame the momenta have following components: 

q = {qo,o,q z ) , 

P = (Po,0,-q z ) , 

I = (lo, l±,lz) > 
^ = (^0> — ^-L' O ' 

k = (k ,k ± ,k z ) , 

k! = (ko, -k±, -k z ) , (3.6) 

where we have used the fact that the quarks are massless. 

Introducing the momentum of the qq pair, K = k + k' we can decompose 
the 3-body phase space of Eq.(3.1) as follows 

dT v dT k dT k ,{2Tif5 i {k + k' + l'-p-l) = -^-dW 2 d$ D d$ 2 , (3.7) 

Z7T 
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where W 2 = K 2 is the invariant mass squared of the qq pair, d<3?2 is the 2- 
body phase space of l+p — > I' + K and d&D is the phase space of K — ► k + A;' 
decay. The explicit expressions read 

d cos uj dx . , 

d * D = ~mT 2^ ' (3 - 8) 

where u and x are the polar and azimuthal angles of k and 

87rA(s,p z , m^J Z7r 

where </? is the azimuthal angle of Z' The function A is defined as A(a, 6, c) = 
>/(a + 6 — c) 2 — 4a6 and s is the total energy squared: 

s = {p + lf. (3.10) 

For given final lepton £' (e~ or u e ) the current matrix element squared 
(without phase space factors) can be written as 

£"'"(A, A'; A B) H^(a, m A, B)D v>li , (A)D^(B) . (3.11) 

A,B= 

7,W,Z 

Here, denote the propagators of electro-weak bosons. The leptonic part 
C flu (X, A'; A, B) is given by the Feynman diagram of Fig. |2] and the hadronic 
one reads 

H^(a, V -,A,B) = e^(p)H^( V -,A,B)e^(p) (3.12) 

with A, B) given by the sum of diagrams shown in Fig.|| for the 

outgoing quark of helicty rj. 

For massless quarks the hadronic tensor satisfies 

ffH^ia, v, A, B) = q u H^(a, rj; A,B)=0 (3.13) 

which allows us to take 

D ^ = WTMi 9 ^- (3 - 14) 

These propagators can be expressed by the polarization vectors by means 
of the following identity: 

et(<l) + £(9) + |4 cS"(?) e£(ff) = ^ - <T ■ (3-15) 
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Using again equations ( |3.13j ) we find that there is no interference between 
different polarizations of the exchanged bosons. Thus we obtain the formula 
for where the contributions from the leptonic and hadronic subprocesses 
are clearly separated: 

A „ _ 1 1 v f dQ^dW 



2 10 vr 4 A( S ,p 2 ,m 2 ) ^ = J (Q 2 + M\){Q 2 + M|) 

7,W,Z 



x Y,L Xp (A,B) H avp (A,B), (3.16) 
p 

where 

L Xp (A,B) = Y,J ^ W(1 ( 9 )e(A,A';AB)e^(g), (3.17) 

H avp (A,B) = J ^d C o S u;e* {p) ^q)H^(a,mA,B)e (p)l/ (q). (3.18) 



The calculation of the leptonic part, L\ p (A,B) is straightforward and 
the result is given by 

L Xp (A,B) = E/^^^{^W^'^)^+^'^)^)^]} 
- E^f^VAp, (3.19) 

A' 

where ei p \ = er p -\(q), g x denotes the coupling of the electro- weak boson B 
to the lepton £ of helicity A, m and m! are the masses of initial and final 
lepton, respectively, and P\ are projection operators: 

P ± = i(i±- 

The explicit expressions for Ly\ p read 



Llr+ = £rl+ = £lr- = ^RL- 





(3.20) 




(3.21) 


2Q 2 Y~(y), 


(3.22) 


2Q 2 Y (y) - 2mm' , 


(3.23) 


—2mm' , 


(3.24) 


2mm' , 


(3.25) 
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where 



Y + (y) = ±, (3.26) 

yZ 

(l-v) 2 

Y-(y) = L ~f L , (3-27) 



2T 

Whl) = 2^ (3.28) 

yZ 



and y is denned as 



y = ¥±£ . (3.29) 

'0 + lz 

In the case of large energy scales we can neglect the electron mass m 
which results in L\/\ p being diagonal in the lepton helicities. Thus for L\ p , 
Eq.(3.19), we obtain 

L Xp (y,Q 2 ;A,B) = ^agfg^L XXp {y,Q 2 ) 

( Y + (y) iovp = X 
= SQ^agfgf { Y-(y) for p = -A . (3.30) 
{ Y (y) forp = 

We have factored out the explicit dependence on a by defining the coupling 
constants scaled by the positron charge gjf = g^/e. 

Let us proceed now with the hadronic part. First we calculate the inte- 
grand of Eq.(3.18). It is given by the Feynman graphs of Fig. ||: 

e\p)M ^ A ' B ) e ( fi )M = G n + G 22 + 2ReG 12 . (3.31) 

Explicit expressions for Gn~ read 

G22 = -/ c ^^%f'^Tr[r^ } (p)/V(p)(9)^W(p)(9)/V(p)(p)], 
Gi 2 = -/ cff ^^'^Tr[^ p) ^ (3-32) 
where f c is the colour factor, g s — the strong coupling constant and 

r = k — p, r'=p — k. (3.33) 

f c is the same for all diagrams and independent of the weak boson type. 
In our case (SU(iV) colour group and fixed final quark flavour) 

fc = \- (3.34) 
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On substituting = (1 + 7775 )/2 into Eq.(3.32) we see that Eq.flinnD 
becomes the sum of two parts — with and without the 75 matrix under the 
traces. In order to perform the integration of Eq.(3.18) we notice that the 
denominators of quark propagators depend only on u and read 

s = (» _ pf = ? - 2„k = _^- + <y + p> + Mw-\eyy oa . 

r- = (tf - P f = / - 2 P y = J* + + p . 

The 'axial' part (containing 75) vanishes upon integration over the angles 
and we obtain the following result for the hadronic part 

H pav = 8/c g 2 s g^ q g* q 'h pa {x, r) , (3.35) 
which depends on two dimensionless kinematic variables x and r: 

PO+Pz 
X = 

qo + q z 

and 



(3.36) 



A(W 2 ,q 2 ,p 2 " 



(3.37) 



W 2 + P 2 + Q 2 

The explicit expressions for h pcT (x,T) are given in Appendix A. In the 
following considerations we will need only their leading terms for P 2 3> Q 2 . 

Before quoting these approximate results let us look at the exact formula 
for the current matrix element squared, Eq.(3.16). Substituting the results 
from Eq.(3.30) and Eq. (|3,35| ) we get 

2VT J A( S ,p 2 ,m 2 ) Zg = J QLn (Q2 + M 2 )(Q2+M 2 ) 

7,W,Z 

x E \ P a p bM 9i q 9v q ' h P°( x > t) > ( 3 - 38 ) 
p y 

where, by definition, 

8iraQ 



n pB M = ^2 L A P (y,Q 2 ;AB)- (3-39) 



r, x, y can be expressed in terms of the integration variables, W 2 and Q 2 . r 
is given by Eq.( 3.37| ) while x and y, defined by Eqs. ( |3.36| ) and ( 3.29j ), read 

A + P 2 -Q 2 -W 2 

X = A + Pi-Qi + W 2 ' (3 - 40) 
2A 

V = 2{s - m 2 ) + A + P 2 - Q 2 -W 2 ' (3 ' 41) 
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where A = A(W 2 ,q 2 ,p 2 ). 

The integration limits defined by kinematics read: 

WL,= (V~s-mf, (3.42) 
Q 2 min = - 2m 2 + ^-(s + m 2 + P 2 )(s + m 2 - VF 2 ) 

max 

=F — A(s,m 2 ,p 2 )A(s,m 2 ,W 2 ). (3.43) 
2s 

At this stage we are ready to discuss the 'equivalent bosons' approxi- 
mation for s > P 2 S> M\,M^,m 2 . To this end we first observe that for 



given P 2 the hadronic part h pa (x,T) is logarithmically dominated by low 
Q 2 values: 



h± T {x,r) ~ ^ 2 log^, (3.44) 

h±±(x,r) ~ (l-x) 2 log^ (3.45) 

with other components finite for Q 2 /P 2 — > 0. Thus we will be interested 
in the region Q 2 <C P 2 corresponding to the leading logarithmic approxi- 
mation. In this limit x gains the meaning of the Bjorken variable of the 
hadronic subprocess: 

p2 

x = — , (3.46) 
2pq 

while y becomes the electroweak boson momentum fraction with respect 
to I: 

pq W 2 + P 2 

y = Vi=^TW- (3 - 47) 

The product of x and y is independent of the integration variables and reads 

* = *y = g. (3-48) 

As shown in Ref.|| the leading terms of h pa are in fact the products of 
the logarithmic factors and the splitting functions P£ (x) and Pq r) (x) of a 
spin-1 boson of helicity p into a quark (antiquark) of helicity rj: 

h p<7 (x,Q 2 ) = ~[P£_» + P|_>)]log^ (3.49) 
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with 

P± ( x ) = p_± ( x ) = 3x 2 , (3.50a) 
P? ± W = Pl(x) = 3(l-x) 2 . (3.50b) 

Changing the integration variable W 2 to y we obtain the formula for 
Ja V which allows for the discussion of the equivalent boson approximation: 



7.W.Z 



X 



" mm 



J/ 2vr 

fK- Q 2 dQ 2 P^ 

Jo* . (Q 2 + M|)(Q2 + M 2) io g Q 2> ^ 



/'/;■ ( - ) • /'/' ' ; 



where 



y max = 1 - 0(m 2 /P 2 ) , (3.52) 
Q mi „ = ^ 2 T ^, (3-53) 
QL x = P 2 Z + V ~ Zy - (3-54) 



P^pBp m * ne aDove equation is a generalization of the electron-boson split- 
ting function, which takes into account the interference effects. 

Let us note that the minimal value of Q 2 is the only place where the 
electron mass must be kept finite to regularize collinear divergencies in the 
case of photon exchange. However, in order to remain within the leading 
logarithmic approximation we have to pay particular attention to the upper 
limit of the Q 2 integration. In general it is a function of P 2 , however 
integration up to the maximum kinematically allowed value Q max would 
violate the condition Q 2 /P 2 <C 1. For our approximation to work we must 
integrate over Q 2 up to some Q max <C P 2 ■ Integrating Eq.(3.51) over Q 2 
within such limits and keeping only leading-logarithmic terms leads to 

J^ = S f - E ^^4B^m»)<»K^,-<r + ^Sr,^C. (3.55) 
A,B,p * 

where F£ B (y, Qmax) is the density matrix of polarized bosons inside the 
lepton and Q 2 is an average Q 2 value. (Exact result from Eq.(3.51) is 

Q 2 = \J OmaxQmin within the leading-log approximation only the order 
of magnitude is relevant). 
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The transverse components of F^ pBp in the case of unpolarized electron 
read (formulae for polarized electron are given in the Appendix B) 

^z + to,« 2 ) - ^tan^ w ^-^ + ' i og g!^, (3.56b) 
fJIz. to, « 2 ) = £ Uo' *w ^ + ^-»»' log , (3.56c) 



,,<> <*. fl pw(l-y) 2 -! , Q 2 + M| , 
F-i+z+iViQ ) = ^tan% — log — , (3.56d) 



^z_(^ 2 ) = £tan% ^-^- y)2 log^S, (3.66e) 



g 1 (i-y) 2 j Q 2 + M 2 
27r4sin 2 e w 2/ ° g 



^ + W + (*,Q 2 ) = ^ZT-ZTTT- log^P^> (3-56f) 



^-w-^) = ^4sk>^*- (3 - 56g) 

where 6w is the Weinberg angle and 

PW= 2^nW" 1 - (3 - 57) 

All other density matrix elements (containing at least one longitudinal bo- 
son) do not contribute at the leading-logarithmic level. 

The density matrix elements F\ B (y,Q 2 ) should be considered as a 
generalized solution to the Eq.(2.9a). At this point we are able to define 
the splitting functions of an electron into a quark at the momentum scale 
Q 2 as 

K~ (Q 2 ) = E tffi 9 E F a p b p (Q 2 ) ® P q P v • (3-58) 

AB P 

This is the generalization of the first sum of Eq.(2.9b). The explicit expres- 
sions for quarks read (see Appendix B for the polarized electron case) 

Pl^Q 2 ) =t{^M*) + *-W]iogg 

n2 i /if 2 

+e 2 tan 4 6>w[$ + (z)+^-(z)]log W + 2 Z 

z 
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O 2 + M 2 1 

-2e 2 tan 2 g w [-^+(^) +Pw^-(^)]log M 2 Z f > ( 3 - 59a ) 

^r^o 3 ) = |f{4[M*)+M*)]kg^ 

+z 2 tan 4 + /4<Mz)] log V ^ 2 Z 

z 

n 2 + _/vf 2 

+2e q z q tan 2 # w [-*-(*) + Pw$+ (*)] log m2 Z 

z 

+ (l + pw) 2 M^)Vlog ^ . ( 3 - 59b ) 



where 



*-(*) = ( 3 - 61 ) 



and 



sm z 6> w 



(3.62) 



with e g and T| being the quark charge and 3-rd weak isospin component, 
respectively. The splitting functions for an antiquark of the opposite helicity 
can be obtained from Eq.(3.59) by interchanging <3? + with <£_. 

These new splitting functions depend on the scale Omax- This scale can 
be in principle fixed by experimental cuts, in such case the final lepton must 
be tagged. A condition commonly used in phenomenological applications 
H is Qmax = eP 2 , where eCl and in general depends on y and zln this 
case the electron splitting functions depend on logP 2 and hence on t. In 
the following we will concentrate on the case with the maximum virtuality 
being independent. In this case the master equations Eq.(2.9b) become 

l F e > {t) = {t) + ^ly y P K (3.63a) 

h'= p=±i 

q,q,G 

h'_= p=±l 
q,q,G 

Note that the convolution of the equivalent boson distributions and 
boson-quark splitting functions, Eq.( 3.58 ), occurs at the level of splitting 
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functions. It is not equivalent to the usually performed convolution of the 
distribution functions because of the P 2 dependence of the boson densities 
Eq.(3.56). Only in the case when the upper limit of integration Q^ax i s kept 
fixed (P 2 -independent), are the convolutions equivalent at both levels. 

To summarize, our splitting functions show two new features. The first 
one, already mentioned before, is the contribution from the interference of 
electroweak bosons (7 and Z only). The second one is their P 2 dependence, 
which results from the upper integration limit Qmax- 



4. Asymptotic solutions to the master equations 

The master equations Eq.(3.63), with spin dependent functions, are most 
conveniently solved in terms of polarized and unpolarized distributions ||: 

4* = Kl+KL (4-la) 



and 



ft = P hl + Pt > (4.2a) 
= Pll ~ Pfc (4.2b) 



= P& + Pll , (4.3a) 



API:' = P£ - P£ . (4.3b) 

Upon substitution of the above definitions into Eq.(3.63) we obtain two sets 
of master equations where unpolarized and polarized functions do not mix: 

J t K" (t) = ±1$ (t) + ^ E Pf ® # (t) , (4.4a) 

h'= 

q,q>G 

I AF? (t) = |- A pf (*) + ^ £ AP "' ® AF ° h > W • (4-4b) 

h'= 
q,q,G 

In the following we will present the asymptotic solutions to these equations 
for the case of unpolarized electron, F£ and Ai^ . 



h+ 

Ji', „ft 
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The equations Eq.(3.63) have been derived within the leading logarith- 
mic approximation which is justified for asymtotically large P 2 values. In 
this region all the logarithms of P 2 over any finite scale are approximately 
equal and in the QCD evolution we choose them to be scaled by Aqcd : 

t = log-^. (4.5) 

QCD 

We present below the solution to the master equations in this asymptotic 
(large t) region. 

The strong coupling constant is approximated by 

a s (t) ~ ^ , (4.6) 

with b = 11/2 — rif/3 for ri{ flavours. Parametrizing the asymptotic solution 
to Eqs.(4.4) for the QCD parton h in the unpolarized electron as 



F e h (z,t) ~ f^{z)t\ (4.7a) 

AFf(z,t) ~ I^ 2 A/f(z)t 2 (4.7b) 

we obtain purely integral equations 

fh = Pf p h ® ffi , (4.8a) 

zo h> 

A/ h as = APf + 1 APH' ® A/ h a f , (4.8b) 
10 h> 

where P% (z) and AP^ {z) are equal to P% (z) and AP^ (z), respectively, 
with all logs put to 1. 

We solve the above equations by a numerical procedure described in 
Ref. |||. We take number of flavours rtf = 5. In the asymptotic region, 
where the quark masses can be neglected, we are left with only two different 
quark distributions: up-type and down-type. 

In figures ||, |5| and || we present the unpolarized quark, anti-quark and 
gluon distributions. The polarized distributions are depicted in figures [7| and 
||. The density zAf^ s (z), not shown here, is more than 10 times smaller than 
zAf^ s (z). One notices significant contribution from the W intermediate 
state in the density of d-type and ii-type quarks. The most surprising 
however is the 7-Z interference contribution which cannot be neglected, as it 
is comparable to the Z term. It violates the standard probabilistic approach 
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where only diagonal terms are taken into account. This also stresses the 
necessity of introducing the concept of electron structure function in which 
all contributions from intermediate bosons are properly summed up. Due 
to the nature of weak couplings they turn out to be nonzero, even in the 
case of gluon distributions. Again the 7-Z interference term is important 
and the W contribution dominates in the asymptotic region. 

One should keep in mind that at finite t the logarithms multiplying the 
photon contribution differ from the remaining ones (cf. Eq.(3.56a)). Being 
scaled by m e , they lead to the photon domination at presently available P 2 . 
The importance of the interference term remains constant relative to the Z 
contribution, as they are both governed by the same logarithm. 



5. Conclusions. 

In the paper we investigated the hadronic content of the electron defin- 
ing a new quantity — the electron structure function. We defined it and 
constructed its evolution equations. We also gave asymptotic solutions to 
these equations demonstrating rich flavour and spin dependence. In all 
cases the parton densities inside the electron grow as logarithm squared 
of the external momentum scale (in our calculation the gluon momentum 
P 2 ). This construction, which looks at first sight a formal manipulation 
only, brings in several new issues. The most important ones are the correct 
treatment of the intermediate electroweak bosons and precise formulation 
of the structure function evolution. In the first case one faces the necessity 
of introducing the electroweak boson density matrix with the off-diagonal 
(7-Z) terms comparable to the diagonal ones. The probabilistic approach, 
one is used to in the partonic picture, can be restored only when looking 
at the whole process i.e. an electron emitting quarks and gluons. Also the 
QCD evolution can be treated correctly only when looking at the complete 
system. The virtuality of the intermediate boson depends in general on the 
external momentum and in this way the electroweak sector enters the QCD 
evolution. The derivation of the electron structure function presented above 
assumes a totally inclusive situation when the outgoing lepton is not tagged. 
In the opposite situation, which is possible only with neutral currents and 
away from the forward direction, the virtuality of the intermediate boson 
can be controlled by experiment and kept independent of the external mo- 
mentum scale. In all cases the known boson structure function should be 
treated with great care in the lepton induced processes. The question be- 
comes very important at momenta where all electroweak bosons come into 
play. But even at presently available momenta, where the 'resolved' photon 
dominates, one can see how the correct treatment of the scales changes the 
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shape of parton distributions. In Fig.|| we compare the asymptotic solu- 
tions of the evolution equations following from our procedure (ESF) with 
those following from naive application of the convolution Eq, (|l.l[) (FF). It 
is possible that the difference can be traced in the analysis of presently 
available data. The inclusion of effects connected with finite momenta and 
construction of the electron structure function which could be applied to 
present experiments is a natural next step along the proposed lines. 
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Appendix A 

Here we present full results for the hadronic part of the QCD current 
matrix element squared as well as some kinematical relations. 
The h pa (x,T) functions of Eq.( p.35| ) read: 



h (x,t) = h ++ (x, t) = (1 — x) 2 C 

16a; 



+ ( * y {-r(l + x) 2 + f (1 - Ux + 49x 2 ) 



{l-xf 



+3f 2 (l - 14x + 25s 2 ) 
+3f 3 (3 - 18s - 19s 2 ) 
+3f 4 (l - s)(3 - 7s) + 3f 5 (1 - x) 2 } C 



(A.l) 



+ y ^ 2 J {8fs(2 - 5s) - 4f 2 (l - 8s + 10x 2 ) 

-6f 3 (l - s)(l - 3s) - 3f 4 (l - a;) 2 } 
-2(1 -x) 2 , 

(s,r) = (x,t)=x 2 C 

+ {-r(l-x) 2 (l+x) 2 

+f(l - 2x 2 - 32x 3 + 49s 4 ) 
+f 2 (l - s)(3 - 13s + 53x 2 - 75s 3 ) 
+f 3 (l-s) 2 (9-38x + 57x 2 ) (A.2) 
+3f 4 (l - x) 3 (3 - 7s) + 3T 5 (1 - x) 4 } C 

+^2 {-8f(l-s)s 2 (l-5s) 

-4f 2 (1 - s) 2 (l - 4x + 10s 2 ) 

-6f 3 (l - s)3(l - 3s) - 3f 4 (l - s) 4 - 16s 4 } , 
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fl - x) 2 

h oo(x,r) = ^ {-r(l + x) 2 + f(l + 2x + 17x 2 ) 

-f 2 (l + 18x - 55x 2 ) + f 3 (5 - 46x + 53x 2 ) 
+3f 4 (l - x)(3 - 7x) + 3f 5 (1 - x) 2 } £ (A.3) 

+ 3(1 2 ~ 2 X)2 {-8fx 2 + 4f 2 x(2 - 3x) 

-2f 3 (l - x)(l - 3x) - f 4 (l - x) 2 )} , 

/i -(x,r) = ho + (x,T) = h-o(x,T) = h +0 (x,T) 

= {r(l - x)(l + x) 2 - f (1 + x - 17x 2 + 31x 3 ) 

+f 2 (l- 23x + 79x 2 - 65x 3 ) 
+f 3 (l-x)(7 - 46x + 55x 2 ) 
+3f 4 (l - x) 2 (3 - 7x) + 3f 5 (1 - x) 3 } C 

+^-rr {+4fx(l - 7x + 8x 2 ) 

-2f 2 (l-x)(l-12x + 19x 2 ) 



-6f 3 (l - x) 2 (l - 3x) - 3f 4 (l - x) 3 + 8x 3 } 



where 



p2 

x 2 Q 2 



T 



(A.4) 



£ = lo gZ2?^' ( A - 5 ) 



(A.6) 



1 -r 

and r and x are defined by Eq.( 3.37j ) and Eq.( 3,36| ), respectively. 
For fixed Q 2 h pa (x,T) are functions of x only with 

x 2 Q 2 

T = P 2 + x 2 Q 2 ' (A " 7) 
Introducing s = s + P 2 — m 2 = 2lp we can define dimensionless variables 

e = ? , (a.8) 

s 

p2 

z = — (A.9) 

s 

in terms of which we can write the formulae needed to change (x, r) variables 
to used in the integral (3.51): 

x 2 f 

r = —\-, (A.10) 

z + x^£ 
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Appendix B 

In this appendix we present the formulae for the boson density matrix 
elements and electron to quark splitting functions for the fully polarized 
case. As can be seen from Eq.(3.51) and Eq. fl3.55 ) the boson density matrix 

elements, F J f B are equal to Pj^g times appropriate logarithms. The 

transverse components of P^ B can be read off Eq.( |3.39 ) and Eq.(3.30) 
and are equal to 

pf pBp = 9±£ e yY± P (y)- (B.i) 

Inserting electroweak couplings we derive following formulae for the trans- 
verse components of B in the case of polarized electron: 



F°Py p (v, Q 2 ) = yY±p(v) log % , (B.2a) 



e 

2 i A/T2 



^~z>Q 2 ) = ^ P 2 wt a n 2 9 w yY- p (y) log Q , (B.2b) 

z 

FlUMQ 2 ) = £ tan 2 e w yY p (y) l og 9l±Ml , ( B .2c) 

z 

< = Z>Q 2 ) = £/>wtan% y y_ p (y) log Q ' , (B.2d) 

z 

<; Z >3 2 ) = -£tan0 wy y p (y) log , (B.2e) 

z 

<w> <J"> = ^-pM log , (B . 2f ) 



F^ pWp (y,Q z ) = 0, (B.2g) 



where 



(l-v) 2 1 
yY_(y)= { - yY+(y) = -. (B.3) 

y y 

The splitting functions of a polarized electron into a polarized quark are 



defined by Eq.( 3.5S| ) with F J f B taken from Eq.(B.2) and a straightforward 
calculation leads to the following result: 

p 9 e =(,,g 2 ) = |M-){^iog^ 

2 2 f An I Q 2 + M| 

+ e 9 p w tan W log — 2 

z 
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£2 , yp. 

- 2e 2 qPw tan 2 % log - ^ Z |> , (B.4a) 

z 



Z7T I y rrig 

+ e^tan 6> w log — — 2 — 
z 



n 2 I 

+ 2e 2 tan 2 W log - ^ 2 Z }> , (B.4b) 

z 



P;i(z,Q*) = ^ + (z) jejlog^ 

2 2 . An , Q'+^Z 

z g/ o w tan (9 w log — — 
z 

Q 2 _|_ 2 

2e q z q p w tan 2 W log — ^ 2 z 

z 



+ (l + m) 2 Vlog^^|, (B.4c) 
^(,,Q 2 ) = |M^){e 9 2 logg 

2. 4 fl , Q 2 + M| 

+ z 9 tan w log — ^ 2 
z 

- 2e g z g tan 2 % log - ^ Z J , (B.4d) 

where $-(z) and p\y are denned by Eq.(3.61) and Eq. (|3,62| ), respectively. 
The splitting functions for an antiquark of the opposite helicity can be 
obtained from Eq.(B.4) by interchanging <£ + with 
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Fig. 1. Lowest order graphs contributing to the process: e + G* — > I' + q + q. 




Fig. 2. Lowest order diagram for D lv ' . 




k 1 k 1 k 



Fig. 3. Lowest order diagrams for 7i v v . 
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Fig. 4. Unpolarized quark distributions zf™ A (z) — solid line. The other lines show 
contributions from different clectroweak bosons. 
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Fig. 5. Unpolarized anti-quark distributions zf^,^(z) — solid line. The other lines 
show contributions from different electrowcak bosons. 
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Fig. 7. Polarized d-quark distribution zAf^(z) — solid line. The other lines show 
contributions from different clcctroweak bosons. 
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Fig. 8. Polarized gluon distribution zA/q s (z) — solid line. The other lines show 
contributions from different clcctrowcak bosons. 
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Fig. 9. Comparison of unpolarizcd d-quark distributions z/| s (z) calculated by 
ESF and FF methods. The upper two lines result from contributions from all 
electroweak bosons while the lower two from 77 only. 



